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1. INTRODUCTION 

Unmanned aerial vehicles are commonly used in a variety of applications, including indoor or 
outdoor surveillance, remote inspection and hostile environment tracking. Quadrotor is emerging as 
a common platform among UAVs owing to its greater capacity for payload and greater maneuverability for 
single-rotor. In addition, quadrotor has a number of benefits in regards to structure, cost and motion control 
compared to traditional helicopters. Quadrotors have therefore received increased attention among 
researchers and engineers, and have also become a promising alternative for multiple unmanned military 
and civilian apps. However, quadrotor UAV has significant science and engineering issues due to its 
characteristics including underactuation, unknown nonlinearities, multivariable and high coupling. Quadrotor 
control is therefore becoming quite complicated and hard, primarily because of its underactuated 
characteristics and nonlinearities. 

It is known that the quadrotor is volatile. Hence, the attitude and altitude stability problems are 
the primary goals for most research in this sector. Some of the methods used to control the quadrotor 
platform are the proportional-integral-derivative (PID) control [1-3], linear quadratic regulator (LQR) 
control [4-6], sliding mode control [7-9], feedback linearization control [10-12], fuzzy logic (FL) 
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control [13-15] and backstepping control [16-18]. In this paper, the quadrotor helicopter's stability issue 
is regarded. The altitude is selected as controllable DOF, along with three angles of attitude, roll-pitch-yaw. 
The dynamic model that describes the movements of the quadrotor helicopter and takes into account 
the different parameters that influence the flying structure dynamics is presented. Thereafter, a backstepping 
approach-based control technique is created. The backstepping control system is a non-linear control method 
based on the theorem of Lyapunov. Due to its recursive design and systematic methodology, 
the backstepping control design methods got many attention [19-23]. Unlike the feedback linearization 
technique with issues such as the precise model requirement, the backstepping strategy provides a choice of 
design instruments to accommodate nonlinearities. The benefit of backstepping compared to other control 
techniques is its flexibility in design because of its recurrent use of Lyapunov functions. The backstepping 
model concept is to pick a suitable variables state recursively to be a virtual inputs for the general system's 
reduced dimension subsystems and for each stable virtual controller, the Lyapunov functions are designed. 
The stability of the overall control system is therefore guaranteed by the designed actual control law. 

Even though the backstepping technique can provide a systematic building process for controller 
design, it is not simple to achieve satisfactory output due to the arbitrary selection of controller parameters 
acquired by the backstepping technique. In order to achieve a good result, it is essential to select appropriate 
parameters because the incorrect choice of parameters can give a poor results and sometime cause system 
unstable. It is also possible to select the parameters appropriately, but it is not feasible to say that the ideal 
parameters are selected. In this work, GSA is used to calculate off-line the ideal parameters for the quadrotor 
system backstepping controller. This study mainly contributes to the backstepping control design approach 
utilizing GSA to maneuver a quadrotor UAV. 


2. QUADROTOR DYNAMIC MODEL 

The quadrotor UAV comprises of four-rotors in cross-configuration, as shown in Figure 1. 
The dynamic equation of movement of the attitude could be deduced from the Euler equation, that could be 
formulated as follows [24]: 
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Figure 1. The quadrotor UAV configuration 
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3. BACKSTEPPING CONTROL OF QUADROTOR 

Only the altitude and angular movement (roll-pitch-yaw) are selected as four controllable 
degrees-of-freedom (DOF) in this study. The dynamic model (1) can therefore be formulated as outlined 
below in a nonlinear dynamic equation: 


X = f(X)+ g(X)u (2) 
where wu the input and X the state vector defined below: 

(u = [u uz ug ul" (3) 

X = [x, x3 x5 x7 |" =[zġ0 4] (4) 


From (1) and (4), the nonlinear functions f (X) and g(X) can be composed appropriately as: 


pa g u,— 0 
Awa, — 8az,Q,q 0 b 0 0 
FO) bas + bade} I=! 4 om o = 
Opas 0 0 bs 


where a, = (ee = tee) [leet a2 = Ir /T ex, a3 = (Izz z lex) /ly > a4 = J;/lyy > a5 = (Lex 7 Tyo) floes 
b, = l/Lx, b2 = l/lyy, b3 = 1/lzz, Uz = (coco ). 

The control objective is to plan a reasonable control law for (2) so that the state direction X can track 
a wanted reference direction X4 = [x1q X3q Xsa X7al’. 

The development of the backstepping control is defined: 

Step 1: Error is defined: 


€i = Xia = %1 (6) 
where x4q is the required path indicated by the reference. The tracking error derivative can then be described as: 
6, = Xyq- %4 (7) 


The first feature of Lyapunov is selected as: 


1 
Vile) = 5 er (8) 
The derivative of V; is: 


V, (e1) = eé, = 61 (žia - *1) (9) 


x, could be seen as a virtual control. The required virtual control value recognized as a stabilizing feature can 
be set as: 


&ı = žia + kiei (10) 


where k, is a +ve gain and ought to be decided by the GSA. 
By substituting its required value for virtual control, (9) then turns into: 


V,(e,) = —k,e? < 0 (11) 
Step 2: The deviation from the required value of the virtual control can be described as: 

e = Žž — Ay = Hy — yg — hye, (12) 
The derivative of e, is expressed as: 


ê, =ï -ú 


= f (x1) + g(x1)uı — Ža — k1ê1 ie 
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The second Lyapunov function is selected as: 


1 1 
V, (e1, €62) = z% + 58 (14) 


The derivative of (14) is: 


V, (e1, e2) = e1, + e262 
= e1 (žia - ž1) + e2% — &1) (15) 
= e, (~ez — k61) + e3 (f (x1) + g(x1)uı — ia — k1ê1) 
= -kef + e,(—61 + f (x1) + ga Ju — ia — kié) 


Step 3: For V> (e1, e2) < 0, the u is chosen: 


1 
u = aa + kii + ia — f (x1) — k2€2) (16) 


where kz is a+ve gain and ought to be also decided by the GSA. Replacing (16) in (15), then: 
V>(e1, e2) = —k,e? = ke? = —E'KE < 0 (17) 


where E = [e,e,]" andK =diag(k,,k2). Since V2(e,,e2) < 0,V2(e1,e2) is negative semi-definite. 
Hence, the system will be stabilized by (16). 


4. OPTIMAL BACKSTEPPING CONTROL 

A controller (16) was intended to stabilize the subsystem in the past section. The coefficients k,, ky 
are parameters of control and must be positive in order to meet stability criteria. These parameters are chosen 
by trial and error in the standard backstepping technique. To overcome this limitation, the GSA is used to 
choose the optimal control parameters value. To select the best parameters of backstepping controller, 
the GSA is used off-line and controller performance differs with adjusted parameters. The quadrotor system 
consists of four subsystems, as mentioned above. Thus, there are eight control parameters to be chosen at 
the same time. An integral absolute error (JAE) is used in this work to assess the controller's performance. 
The index JAE is defined as below [25]: 


t 


IAE = flew dt (18) 


0 


Since the system is consists of four subsystems, a vector JAE for the whole scheme is therefore considered 
as IAE, = [IAE, IAEp IAEp IAEy], where the subscripts are signified for altitude, roll, pitch and yaw 
subsystem, respectively. GSA aims at minimizing the fitness function J, stated as: 


J = IAE;-W (19) 


Where W=[W, W, W; W,]’ is weighting-vector utilized to establish the importance of 
the several parameters and “W” changes in the range of 0-1. Same weights for the four objectives to be 
achieved are considered in our case to indicate the error indexes minimizations are correspondingly essential. 
The time-domain simulation is performed for the simulation period, t for the fitness function calculation. 
This fitness function needs to be minimized in order to expand the system response in terms of the settling 
time, overshoots and steady-state errors. 


5. SIMULATION RESULTS 

The performance of the suggested strategy is assessed in this section. The parameters values of 
the quadrotor used in the simulations are given as [26]: m=0.5 kg, /=0.2 m, g=9.81 m/s’, [,=Iy=4.85x10" 
kgm’, [,,=8.81x10° kgm’, b=2.92xx10° Ns’, d=1.12xx107 Nms”. Then, GSA is utilized to tune the backstepping 
control parameters. The following values are allocated in this research for optimization of controller parameters: 
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The number of agents (masses)=15; 

The search space dimension=8 (i.e., Kj-1.g); 

The parameters searching ranges are restricted to [0, 15]; 
The number of maximum iteration=20; 

Optimization process is repeated for 20 times; 

The simulation time, t is equal to 10s; 

The best optimized controller value is selected among the obtained finest value. Figure 2 shows 
the range of the fitness function and the iterations number. In the meantime, the evolution of the control 
parameters is shown in Figure 3 over the iteration number. The GSA technique can prompt convergence 
and achieve desired fitness value through about 20 iterations. 
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Figure 2. Fitness function with iterations number 
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Figure 3. Control parameters versus iterations number 
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To investigate the efficacy of the suggested optimal backstepping controller, two simulation studies 
were conducted on the quadrotor. The first simulation is conducted for stabilizing problem and the second 
is performed to investigate attitude tracking problem. 


5.1. Simulation experiment 1: stabilizing problem 

The altitude/attitude control experiment is implemented in this simulation. The quadrotor is required 
to achieve and hold at a desired altitude/attitude. The desired altitude/attitude are set 
at Xiq = [Za Pa Oa Pal = [15,0,0,0]". The initial states are set to z = 0, @ = 0.2,0 = 0.2 and w = 0.2. 
By using the proposed OBC, it can be seen that the quadrotor can fly toward the desired altitude and manage 
to maintain at the desired altitude/attitude as shown in Figure 4. It can be also observed in Figure 4, 
the attitude converges quickly to stabilize the quadrotor. 
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Figure 4. Altitude/attitude response 


5.2. Simulation experiment 2: attitude tracking problem 

The attitude tracking control experiment is conducted in this simulation. The quadrotor is required 
to track a desired reference signal. The sinusoidal signals are utilized as a trajectory to the attitude angles. 
It can be observed in Figure 5, the quadrotor manages to follow the trajectories efficiently. The findings also 
indicate that, the OBC manage to give a satisfactory tracking performance since the tracking error is small. 
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Figure 5. Attitude tracking response 
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6. CONCLUSION 

In this paper, the application of an optimal backstepping controller for maneuvering a quadrotor 
UAV is successfully demonstrated. First, the quadrotor's mathematical model is introduced. The optimal 
backstepping controller is then developed in which the controller can select the parameters automatically via 
GSA. To ensure the system's stability, the backstepping control design is based on the Lyapunov function. 
Finally, the OBC is applied for a quadrotor UAV stabilization and trajectory tracking missions. The results 
indicate that the suggested control scheme can achieve a good transient and tracking response. 
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